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Abstract 

We reconsider the issue of superluminal propagation in the DGP model of infrared mod- 
ified gravity. Superluminality was argued to exist in certain otherwise physical backgrounds 
by using a particular, physically relevant scaling limit of the theory. 

In this paper, we exhibit explicit five-dimensional solutions of the full theory that are 
stable against small fluctuations and that indeed support superluminal excitations. The 
scaling limit is neither needed nor invoked in deriving the solutions or in the analysis of its 
small fluctuations. To be certain that the superluminality found here is physical, we analyze 
the retarded Green's function of the scalar excitations, finding that it is causal and stable, but 
has support on a widened light-cone. We propose to use absence of superluminal propagation 
as a method to constrain the parameters of the DGP model. As a first application of the 
method, we find that whenever the 4D energy density is a pure cosmological constant and a 
hierarchy of scales exists between the 4D and 5D Planck masses, superluminal propagation 
unavoidably occurs. 



Introduction 

The possibility of modifying general relativity at large distances is both theoretically fascinating 
and of potentially great import for cosmology. The theoretical challenge stems from decades of 
negative results about the very possibility of changing the long distance behavior of gravity. 
Various roadblocks stand in the way of such a modification. The first is the celebrated van 
Dam-Veltman-Zakharov discontinuity [I], which is due to the fact that in flat space a massive 
spin-2 field propagates a zero helicity state. In the massless limit, that helicity becomes a scalar 
that couples with gravitational strength to the trace of the stress-energy tensor of matter. Such 
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coupling changes the GR prediction of, say, the bending of light by massive objects by a finite 
amount incompatible with experiment. This is to say that the massless limit of a theory of massive 
gravity is not GR but rather a Brans-Dicke theory with Brans-Dicke parameter uj = 0(1). Current 
solar-system observations imply uj > 40,000. More theoretical but potentially worse problems are 
the existence of ghosts around Schwarzschild-like solutions [3], the absence of a lower bound 
on the energy of small fluctuations around generic backgrounds [I], and the existence of an extra 
propagating degree of freedom, beyond linear order, in generic nontrivial backgrounds (also in jl]). 
All these problems are related, and call for a more sophisticated way of changing gravity at large 
distances. The DGP model [5] provides at least in part a cure for some of the problems just 
mentioned. In DGP, gravity is four dimensional below a certain distance L, and five dimensional 
above L. This distance is given by the ratio of the 5D Newton's constant over the 4D Newton's 
constant, L = G^/G^^ Thus, the familiar r~ 2 behavior of Newton's law becomes a weaker r~ 3 
at large distances; empirically, L = C(-RHubbie)- While still suffering [U] from the severe strong- 
coupling problem which also plagues massive gravity [7] at the quantum level, the DGP model 
does not propagate ghosts on asymptotically flat backgrounds [8]. Ref. [8] also shows that the 
crucial aspect of the DGP dynamics, namely the decoupling of the scalar degree of freedom around 
a massive source, holds in a decoupling limit in which the 4D Planck mass M4 — > 00, while the 
quantum DGP scale A 3 = (M4/L 2 ) 1 / 3 is kept constant. In this limit, gravity around a source 
of (fixed) Schwarzschild radius Rs is in the linear regime whenever r 3> Rs, with the exception 
of a self-interacting scalar mode, 7T, that couples with gravitational strength to the trace of the 
stress-energy tensor. The strong self-coupling of the scalar causes it to decouple at distances r <C 
(RsL 2 ) 1 / 3 P, [101 HI [8], and to propagate superluminally in the range R s < r < (i?sL 2 ) 1/3 P2]- 
A recent work [13] shows that this pathology of superluminality is generic in the decoupling limit. 

A subtle aspect of the decoupling limit is that it washes away the non-local part of the scalar 
propagator. In DGP, the degree of freedom that becomes the scalar of the decoupling limit is one 
of the helicities of the graviton, whose propagator is non-local. The DGP graviton is indeed a 
resonance, not a stable particle (it can decay from the 4D brane into the 5D bulk). Equivalently, 
its propagator (better, its retarded Green's function) does not have poles on the physical sheet 
of the complex frequency plane, but only cuts along the real axis. Moreover, the kinematic 
region of interest for studying locality is precisely that where the non-local part of the propagator 
cannot be neglected. Another possible complaint about the conclusions reached in [121 US] is 
that the background on which the scalar mode propagates is only a solution of the decoupled 
equations and not of the full DGP model. Finally, in a model such as DGP in which 4D non- 
localities stemming from the integration of the 5D bulk physics exist, one must exercise care 
in distinguishing true superluminal propagation from non-localities introduced by an incorrect 
definition of the sources [i~4"j . 

To respond to all these doubts, after a lightning review of the DGP equations in section 2, 
we will review some of its exact solutions in section 3; the time independent ones of [6] and the 
time-dependent ones of [T7j . In particular, we compute the kinetic term of the scalar DGP mode, 
7r, in the backgrounds found earlier. From that computation we find a speed of propagation of 
signals cv greater than c, thereby confirming the results of [Bj. Section 4 is devoted to studying 

1 Equivalently L = M%/M$, where M4 is the 4D Planck mass while M5 is the Planck mass in 5D. 
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the issue of non-locality of rr and of a proper definition of its sources. To be sure that it indeed 
propagates superluminally, we study its retarded Green's function without resorting to the local 
approximation. Our analysis confirms that the Green's function is indeed nonzero inside an 
enlarged light cone, determined by c^. Equally importantly, we find that the Green's function 
does not have instabilities growing exponentially in time, in accordance with the analysis done in 
the decoupling (or local) limit. 

In section 3 we tackle a di?erent problem. We examine whether a range of parameters exists, 
for which the DGP model admits only solutions without superluminal propagation of signals. In 
other words, we try to use special-relativity causality to constrain the parameter space of DGP. 
We find that, at least in the case that the 4D brane energy is a pure cosmo logical constant, no 
range of parameter is safe. This may point to an upper bound on L (the DGP length scale) or on 
the horizon radius in any relativistically invariant, causal completion of the model, such as those 
provided by D-brane or other constructions in string theory. 



1 The model, and our goal 

DGP is a 5D model where gravity propagates throughout an infinite 5D bulk, and matter fields 
are confined to a 4D boundary. The action for gravity at lowest order in the derivative expansion 
is a bulk Einstein-Hilbert term and a boundary one, generically with two different Planck masses 
M5, M4 (plus a suitable Gibbons-Hawking term). Even if we set the gravitational boundary action 
to zero at tree-level, matter loops will generically generate a localized 4D Einstein-Hilbert term. 
The phenomenologically interesting regime is M 5 <C M 4 . 

The DGP equations are the Einstein equations in the bulk, q 

G^mn = J (1) 

along with the Israel junction conditions on the brane (with Z 2 symmetry across the brane), 

M 2 4 G, U - 2Ml (K^ - g^K) = (2) 

(a boundary is equivalent to a brane with Z 2 symmetry). is computed with the normal vector 
pointing into the bulk. Note that we still have V M T Ml/ = 0, by virtue of (K^, — g^K) = 0, 
coming from the momentum constraint implicit in Q. 

Ref. [6] studied small fluctuations about a generic curved solution <]mn- Since the stress-energy 
tensor is localized on the boundary, it is interesting to 'integrate out' the bulk modes and rewrite 
the theory as a 4D theory. Namely, one can solve the linearized 5D Einstein equations for given 
sources and given 4D boundary conditions, plug the solution back into the action, and thus get an 
effective four dimensional action that only involves the sources and the 4D field configuration on the 
boundary. The downside of doing so is that such an effective action is non-local, because massless 
degrees of freedom have been integrated out. Remarkably, at short-distances the effective action 
is approximately local — hence the statement that in the model gravity looks four-dimensional at 

2 We are using the mostly plus signature for the metric. Also we define the D-dimensional Planck mass as 
Mg~ 2 = l/8nG D . 
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short distances. Indeed the resulting quadratic action for small fluctuations at distances much 
shorter than the solution's curvature radius and than the critical DGP length-scale is [6] 

~ Jd A x^g\Ml{\h^Uh^-\hUh) + \h^5T^ 

+ 3m 2 M 2 tt Qtt + 2mM\ tt (K^ - g^K) V M V u tt +\mix 5T^ J , (3) 

where m = 1/Ldgp — 2M^/M% is the DGP mass scale, all 'barred' quantities are computed on 
the solution g^ u , and we are implicitly assuming de Donder gauge for h^ u : V M (/i Mi , — \g^u h) = 0. 
7r is a scalar degree of freedom, which can be thought of as a 'brane bending' mode. Indeed in one 
gauge it measures fluctuations in the position of the brane along the fifth dimension. In another 
gauge it is a component of the five-dimensional metric. The gauge-invariant, physical statement 
is that it couples directly to T^, and thus mediates a physically measurable forcej^] The non-local 
contributions alluded to above, are corrections to the above action schematically of the form 



MimhV-Dh , M 1 ^ tt \J -Utt . (4) 

At high momenta, d ^> m, their effect is negligible with respect to the local kinetic terms we 
explicitly kept. For the moment, we will concentrate on the local part of the quadratic action, and 
come back to the interpretation (and a more careful definition) of these non- localities in Sect. [3j 

We see that the tt kinetic action — thus its propagator — depends on the extrinsic curvature 
of the solution. In the following, we will look for 5D background solutions where this extrinsic 
curvature correction makes tt excitations superluminal while keeping their energy positive. That 
is, we will look for stable solutions with superluminal excitations. 

An important remark is in order. One might think that, even if we find superluminal excitations 
of curved solutions, these excitations will simply correspond to signals connecting two points A 
and B on the curved boundary by taking a straight shortcut through the bulk. In this case our 
4D superluminality would be fake — no signal would exit the true 5D light-cone. This however 
cannot be the case, for two reasons. First, the extrinsic curvature correction to 7r's kinetic action 
above is local: it can give a finite shift in the local propagation velocity of tt, for arbitrarily close 
A and B. The shortcut effect, instead, is only efficient for A and B far from each other, so that 
the curvature of the boundary becomes visible. The shortcut fractional correction to the average 
velocity goes to zero if we bring A and B closer and closer. In other words, a shortcut effect 
would correspond to a non-local term in the effective 4D Lagrangian. Second, linear tt excitations 
correspond to pure gauge deformations of the bulk geometry [6]. Therefore, when we excite tt no 
physical information is traveling through the bulk. We can choose a gauge such that the bulk is 
unperturbed, and the geometric excitation is manifestly localized on the boundary. 



2 Stable solutions with superluminal tt 

In order to (a) find bulk solutions, and (b) study the kinetic action for tt about them, it 
is convenient to focus on configurations with a high degree of symmetry. One possibility is to 

3 To get to ([3]), tt has been demixed from the metric perturbation h^ v j6]. This is the origin of the coupling 
7rT M M . Equivalently, matter couples universally to the 'Jordan-frame' metric perturbation = + rair g^. 
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consider spatially flat cosmological solutions, for which the bulk geometry is exactly flat. There, 
however, one can show that imposing the null energy condition (NEC) for T^ u and stability under 
7r fluctuations, automatically implies stzMuminality of ix [TB]. The next-to-simplest possibility 
is solutions that are spherically symmetric (in five dimensions), with the brane curled up as a 
3-sphere. In such a case, Birkhoff's theorem forces the bulk metric to be static, and of the 5D 
Schwarzschild form: 

*fa = -f» dt2 + J^dr 2 + r 2 dn 2 3 , f(r) = (5) 

(here and in the following we use f 2 to denote what is usually called /). The parameter p is 
related to the solution's 5D ADM mass by 

M ADM = 37r 2 M 5 V . (6) 



2.1 Time independent spherical solutions 

We first look for solutions with a static spherical brane located at some fixed value of r. This 
case has been considered in [6]. From the 4D viewpoint, the brane is a static Einstein Universe. 
The 4D induced metric is 

g^dafdx" = -f{r)dt 2 + r 2 dtt\ , (7) 

and the extrinsic curvature is ^ 

K^ = -fd r g^. (8) 

We have chosen the bulk to lie outside the brane, so that it is asymptotically flat and free of 
horizons and singularities. Then the DGP equations (|2| read [6] 



1 

(mr) 2 

where m = 2M^/Ml is the DGP mass parameter, r is the brane position, / is computed at r, and 
p and p are the 4D energy density and pressure. Note that the energy-momentum conservation 
equation is trivial in this time-independent case, so we must use both the space-space and time- 
time components of the Einstein equations. 

At short distances the quadratic ix action (j3| is 

C^M 2 [3m 2 nDn + 2mrr(K^-g^K)V fi V u n] , (11) 

so the new kinetic matrix is proportional to 

3m 2 9tw + 2m{K„ v - g, u K) = 3m 2 diag (- [l - ^] f 2 ; [l - ^(/ + 1//)] 9ij ) . (12) 



1 
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M 2 m 2P ' 



(9) 
(10) 
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If 7r is to have positive kinetic and gradient energies, we want both terms in brackets to be positive. 
For 7r to be subluminal, we want the second bracket to be smaller than the first. Indeed, the speed 
of 7T excitations as measured in a local Minkowskian frame at rest with respect to the brane is 



C 7T — 1-2/ ' ^ ' 



as can be seen by locally going to Riemann normal coordinates & 

dp = fdt 2 , df 2 = g ij dx i dx j = r 2 dtt\ . (14) 

Specialize now to the case of a pure cosmological constant on the brane, p = —p. From 
equating Q and (10) we find 

" + — (\ ~ 2f) = 0. (15) 



(mr) 2 mr \f 

Given any value mr > 0, this determines /, thus m 2 /i. The cosmological constant M f 2 is then 
determined through either (|9| or (10). Thus, there is a single parameter family of possible static 



solutions, one for each r > 0. 

The 5D mass p is positive only for mr > 2. We want p to be positive, otherwise the solution's 
5D mass would be negative, thus signaling an instability against accelerating away. The cosmo- 
logical constant is always negative, ranging monotonically from — 3/2 m 2 M| to as r ranges from 
to oo. Notice that the brane always lies outside the would-be Schwarzschild horizon. 

Figure [T] shows the stability and superluminality of the solutions. The n mode has ghostlike 
(negative kinetic energy) and/or tachyonic (negative gradient energy) instabilities for mr < 2, 
corresponding to p < 0. It becomes stable thereafter, but is always superluminal. This means 
that in DGP, for values of the brane cosmological constant in the range -3/2 m 2 M| < p < 0, there 
exists a solution which can not be thrown out on stability grounds, which nevertheless supports 
superluminal propagation. As r — ► oo, the kinetic and gradient coefficients approach 1 from below, 
and the speed of 7r approaches 1 from above. Figure 2 plots the 4D energy density p vs. the brane 
position r. 



2.2 Time dependent spherical solutions 

By studying static spherical solutions, we have seen that for a finite range of negative values of 
the brane cosmological constant, — 3/2 m 2 M| < p < 0, there exist stable superluminal solutions. 
To obtain information about other values of the brane cosmological constant, we now search for 
brane solutions that are still spherical, but with a time dependent radius r = a(t). For a 4D 
observer, these solutions describe a cosmology with positive spatial curvature. These backgrounds 
have been studied with somewhat different methods first in [15] and subsequently in [17] . Here 
we will study the stability and speed of propagation of their fluctuations. 

To start, take coordinates t,9i as brane coordinates. 6i, i = 1 . . . 3, are the angles for the 
three-sphere. The embedding is then 

t = t, r = a(t), 9i = Bi. (16) 
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Figure 1: Stability and superluminality of static solutions. The thick black curve is the one-parameter 
family of solutions discussed in the text. The upper-left white region is where the brane is inside the 
Schwarzschild horizon — and of course there is no static solution there. 



The induced metric is 



ds 2 = - 



P 



da 



dt 2 + a(t) 2 dQ 2 3 . 



(17) 



Note that the term in brackets in front of dt 2 should be positive, to ensure that the brane is 
time-like. This requires |^jf | < f 2 , i-e. that the brane not be expanding or contracting as fast or 
faster than the speed of light in the bulk. 

The tangent vectors to the brane are e A [1 = 7^-, where X A are the 5D embedding coordinates, 
and y a are the brane coordiantes. We have, 



da 

1,^,0,0,0), 



(0,0, e,). 



(18) 

(19) 
(20) 



The normal vector n A is orthogonal to all of these, is spacelike, and is normalized, i.e. qab 
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Figure 2: p as a function of r for static solutions. 



0, g A BU A n B = 1. 



n 



I / da 4 



da \ 2 , m \ (it 



Kr. A 0,0,0 . 



/v-(f) +/ 



Again, the bulk is taken to be outside. 



(21) 



The extrinsic curvature can be calculated using K^ v = e A ^e B v V a^b = e. A d v nA — eA u eB v^AB n c, 



K, 



fJ,U 



m 2 +f 4 



dn 2 {3) . (22) 



Here f = £ . 

We want to bring the induced metric into FRW form. To do this, define a new time coordinate 
r(t) on the brane by solving 



/ dr\ 2 2 I f da 
\dt) = * ~J*\dt 



The embedding is now 
and the induced metric is 



dr f 

t = t( r ) r = o(r) , 0* = ft 



,/f i a „ , / a +r$v. 



\dr J 



ds 2 = -dr 2 + a(r) 2 dQ 



3 ' 



(23) 

(24) 
(25) 
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Some quantities we will need are 



^00 

Gin 



R'ij g R 9ij 



a 2 a 2 



- (l + a 2 + 2aa) 7^ 
(1 + a 2 ) 7fj . 



(26) 
(27) 
(28) 



Dot means -4-, and 7^ is the metric on the unit three-sphere. The spatial part of the extrinsic 
curvature, in these coordinates, reads 



Kij = a^f 2 + a 2 



la 



(29) 



(There's actually no need to calculate Kqq, we can get away without it, as we'll see.) 
Trace reversing pj), we have 



Mi I R„ 



2MSK 



1. 



V*> g ± v if^ 

Using the quantities above, the spatial components give us the Friedmann equation, 

a 2 1 \ m — 1 

a 



a 2 a 2 



3Mf ' 



where 



m 



2M| 



(30) 



(31) 



(32) 



is the DGP scale and p the total 4D energy density. We also have the usual conservation equation 



p + 3-(p + p) = 



(33) 



which together with (31) determines the 00-component of (30). (Note that in the static case, the 
conservation equation is trivial, so there we did need the 00-component of (30) as well.) 
We can find the on-shell value of K 00 by looking at the trace of (|2]), 



M 2 R + 6MlK = T 



-p + 3p 



which implies 



K 



m 
3 



3M| 



fa 2 


1 d\ 






\a 2 


a z a J 



a 



The coefficient of the 7r kinetic term reads 

,2 



3m + 2m(K llv - g^K). 



(34) 

(35) 
(36) 

(37) 



9 



n 2 



The OO-component is 

-3m 2 \\-—y/p^ 
ma 

The bracket must be positive if there are to be no ghosts. The r/'-components are 

2 



3m 2 a 2 'jij 



3ma 



(38) 



(39) 



The bracket must be positive if there are to be no tachyons. The local speed of tt is 



1 - 



2 

3ma 



a 2 



2 

3m? 



3Mj 



+ 2! 



(40) 



Specialize now to the case of a pure cosmological constant on the brane, — p = p. Differentiating 



the Friedmann equation (31), we find an equation which can be solved for a in terms of lower 
derivatives. There is a three parameter family of solutions; we can pick arbitrary initial data 



a(0) and d(0), and some value for p. The cosmological constant is then determined by (31). The 



time derivative of (31) is used to evolve these initial data, creating a trajectory in the space of 
independent initial data spanned by a, a and p. Note that the initial data must also satisfy 
1 - 



+ a > because of the square root in (31 ). 



Figure [3] shows the a, a, p parameter space. The plot shows the region where tt is stable and 
superluminal. We will try to rule out more values of the brane cosmological constant by asking 
whether they support solutions which stay in this region. In addition, the solutions must be within 
the effective field theory where tt non-linearities can be ignored^ 

We seek solutions with positive 5D energy so we restrict our attention to the region p > 0. 
The value of p is fixed and conserved, so all the dynamics occurs within a p — const, plane. This 
plane is shown in figure |4j Also shown is the region where tt is stable and superluminal. We seek 
solutions that stay in this region for all times. 



The Friedmann equation (31) gives another surface with a fixed value of p on which the 
dynamics is confined. The intersection of these surface with the p = const, plane then determines 
the dynamical path. The time parametrization along the path is determined through the evolution 
equation given by differentiating the Friedmann equation with respect to time and solving for a. 
The evolution can in fact come to an end at points where a = oo. This can happen at points 
where a, a are finite, so this represents a 4D curvature singularity, as can be seen from the Ricci 

scalar R — 6 f| + ^ + (these are the pressure singularities discussed in [IT]). As it turns out, 

the line on which a blows up is exactly the border of the stable superluminal region. This means 
that solutions can never cross into or out of this region. 

For any value of p there are solutions which live out their lives in the stable superluminal 

(l-y/T+- 



region, ror p > — - 



2/i 



m 2 p 



they come into existence out of a pressure singularity at 



incidentally, the region in which tt is both stable and subluminal does not intersect the region [i > 0. This 
means that in this class of solutions there are none which, at any time in their history, are subluminal, stable, and 



have positive energy from the 5D point of view (^ > 0). 
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one border of the stable superluminal region, travel through the region, and then come to an 

3M 2 ( I \ 

end in a pressure singularity on the other side of the region. For p < \ 1 — y 1 + m 2 p), the 

solutions have a pressure singularity in either the past or the future, and live for an infinite proper 
time, either coming in from or going out to a = oo. Thus, for any value of p, there exist solutions 
which cannot be thrown out on stability grounds, and yet have superluminal propagation. 

As a consistency check of our results, we want to recover the celebrated 'self-accelerating' 
solution — a zero-tension de Sitter brane embedded in a Minkowski bulk. Therefore, we look for 



an empty brane (p = 0) in a flat bulk (p = 0). The Friedmann equation (31) reduces to 



1 + a 2 = ma Vl + a 2 . (41) 
Integrating the Friedmann equation yields the unique (up to time translations) solution 

air) = — cosh(mr). (42) 
m 



Rewriting in terms of the coordinate time t, using (23), we obtain 

1 



a(t) = —Vl + mH 2 . (43) 
m 



u-d Schwarzschild horizon 




Figure 3: a, d, p parameter space. 



11 



This is the familiar self-accelerating solution, a hyperboloid in a flat bulk, where the bulk is taken 
to be the outside. The hyperboloid is maximally symmetric, embedded in flat Minkowski space, 
so its extrinsic curvature is proportional to the induced metric, 



K. 



and the scalar kinetic matrix (37) becomes 



3m 2 g^ u + 2m(K flu - g^K) = -3m 2 g^ u . 



(44) 



(45) 



5-d Schwartzschild horizon p ~>p c 

i I 



a 



-1 



-2 



0.0 



0.5 



1.0 




a = cc 



Figure 4: Examples of stable, superluminal solutions. Shown is a typical p = const. > plane. The 
solution is confined to move along a p = const, contour, shown here as solid black lines for two different 
values of p, with arrows indicating the direction of flow. The solutions come to a 4D curvature singularity 
when a = oo, which is exactly the boundary of the region in which it is stable and superluminal. There is 
a critical value p c = fl — \A + m 2 p S j , at which the p = const, line bifurcates into two pieces. The 
leftmost line shown is for a value p > p c , and the stable superluminal solution lives for only a finite proper 
time, going from singularity to singularity as it traverses the stable superluminal region. The rightmost 
line is at a value p < p c , and the stable superluminal solution lives for an infinite time, extending to 
infinite radius. The values used for this plot are, in units of m, u = 1.5, p = ±0.6M|. 
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Thus the extrinsic curvature corrections have the effect of reversing the sign of the kinetic term, 
recovering the well known fact that tt is a ghost, propagating at exactly the speed of light [6]. 



3 The retarded Green's function 

We now come back to the issue of non-local corrections to the 4D effective action ([3]). 
In particular, we want to determine the support of 7r's retarded boundary-to-boundary Green's 
function. This is enough to study the causal response of n to local boundary sources, i.e. to study 
the causality properties of the interaction between sources mediated by tt. We will see that to this 
purpose, the conclusions based on the local part of the action are correct. 

To simplify the analysis, consider a DGP-like theory for a scalar, with a localized superluminal 
(c > 1) kinetic term on the boundary: 

□ 5 tt(x, y) + 5(y) □ 4 vr(x, y) = J{x)5{y) , (46) 

where 

□ 4 = -d 2 + c 2 V 2 (47) 

and J(x) is the boundary source. For simplicity we are suppressing the dimensionful constant m 
that sets the critical scale. So, here the critical scale is one. This example is much simpler than 
DGP, but it captures the features that are relevant for our discussion. 
Taking the Fourier transform in 4D, the Green's function 

G(x,y) = J ^i^G(p,y) (48) 

obeys the equation 

[(p 2 -d*) + 5(y)p 2 ]G(p,y) = 5(y), (49) 

with 

p 2 ee -to 2 + c 2 p 2 . (50) 
The solution that vanishes at y — > ±oo for Euclidean four-momenta, p 2 > u 2 , is 

g(p, y) = — 9 r^= 2 ■ (51) 

—uj z + c z p z + lyp — oo z 

In the complex u plane, for real p this function has branch cuts starting at u = ±|p|. We define 
the square root in such a way that the branch cuts run along the real axis from the branch points 
to u — ±oo, respectively. This corresponds to setting the branch cut of \fz along the negative 
real z- axis: 

Vpe 10 ee ^/pe ia/2 , for - tt < a < tt (52) 
G(p, y) has poles too, for 

- uj 2 + c 2 p 2 + 2y/p 2 -u 2 = . (53) 
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Figure 5: The integration contours for the integral in uj, as described in the text. The branch points are 
at u = ±k = ±|p|. 



However it is easy to convince oneself that in the superluminal case c 2 > 1 these lie on the second 
sheet — i.e., they correspond to zeroes of the denominator in (51) for the other branch of the 
square-root. To see this, define 

z = \/p 2 - ^ 2 • (54) 
Given our definition of the square root, z takes values in the right half-plane, Ke(z) > 0. Now, 



the pole condition ( 53 ) reads 

z 2 + 2z + (c 2 - l)p 2 = , (55) 

which for real p and c 2 larger than one, has no solutions in the right half-plane. 

The fact that there are no poles in the physical sheet of u tells us that there are no free so- 
lutions for 7r with the given boundary conditions at \y\ — > oo. This is hardly surprising, for we 
already know that in DGP there is no mode localized on the boundary. Rather, four-dimensional 
excitations behave like resonances, "leaking" into the bulk with a timescale of order m~ x . Never- 
theless, the Green's function for it is all we need to compute the interaction between local sources 
that couple to it, regardless of whether it describes an asymptotic state of the theory. 

The Green's function in real space is 



G(x,y) 



d A p 



-iujt+ip-x— y/p 2 — ^ 2 \y\ 



(56) 



-on 2 + c 2 p 2 + 2\Jp 2 — uj 2 

To define the retarded Green's function, we impose the condition that G vanish for t < 0. If we 
integrate just above the cuts on the real u axis, as shown in fig. |5]a, we can close the contour at 
infinity in the upper half-plane, where we have no singularity. At infinity the integrand behaves 
like 

e-**e-*M , (57) 

where, as we discussed above, z has a positive real part. This means that G(x,y) vanishes for 
t < 0, as desired. In conclusion, the retarded Green's function corresponds to integrating above 
the cuts in the u plane. 
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We want to see what happens for t > 0. First, notice that the retarded Green's function 
vanishes in the bulk for \y\ > t. This has to be expected, for we do not anticipate any superluminal 
effect in the direction orthogonal to the boundary. To see that this is indeed the case, suppose 
we still compute the cj-integral by closing the contour in the upper half-plane. The integrand's 
behavior is still given by (57), with t > 0. For fixed real p and large u in the upper half plane we 
have 



\uj\ > \p\ , Im(a>) > 



z — > — iuj 

e ~iut e -z\y\ 



^l-p 2 /uj 2 
D iw(|y[-t) 



(58) 
(59) 

(60) 



which decays as long as \y\ > t, and we get 

G Tet (x, \y\>t) = , 
as expected. 

Notice that in the case of tt, what happens in the bulk is not physically relevant, since 7r is pure 
gauge there. More interesting, and more complicated, is the propagation parallel to the boundary. 
We focus on the y = case, which corresponds to the boundary-to-boundary propagator. First, 
we can integrate over the angular variables of p. Defining k = \p\ and r = \x\ we get 



G mt (x,y = 0) 



1 1 

27r 3 r 



dk k sm(kr) / duj 



-iut 



-uj 2 + c 2 k 2 + 2 V / F^" 



(61) 



Then, for t > 0, we can integrate in uj by deforming the contour in the lower half-plane as shown 
in fig. ^p. The contribution from infinity vanishes, and we are left with the integral along the cuts: 

1 1 f°° 

dk k sm(kr) x 



G ret (x,y = 0) 



2vr 3 r 
-fc 



du 



duj le-^'A 



1 



-u 



2 + c 2 k 2 + 2VW 



(62) 



where A[- • •] denotes the discontinuity across the cut — the value above minus the value below. 
For our definition of the square root, we have 



k < uj < oo = 


> Vk 2 -uj 2 




\J k 2 — uj 2 


oo < uj < —k = 


> Vk 2 -uj 2 



Also 



e = uj 2 - e 



vJ- 



above 
below 
above 
below 



-K 

+i( 
-K 



with ( G 
with C G 



(d( = uiduj . 



(63) 
(64) 
(65) 
(66) 

(67) 



The two uj- integrals in (62) thus combine to yield an integral all along the real £ axis, 

„—lLjt 



G iet (x,y = 0) 



1 1 

2ir 3 r 



dkksin(kr) / d( 



o 



2 1 

7r 3 r 



dkksin(kr) / d( 



sin(a;t) 



u [C 2 - (c 2 - l)k 2 } +4C 2 

c 2 



u [C 2 - (c 2 - 1)A; 2 ] +4C 2 



(6J 
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where for real £ 

cu = sign(C) v / C 2 + A; 2 , (69) 

and in the second step we only kept the even part of the integrand under ( — > — £. What is 
interesting about this form of G Tet is that sm(ut)/u>, as an even function of its argument, is 
analytic in the whole ( and k planes, and we do not have to worry about the cuts of uj as a 
function of ( and of k. 

In the k plane, we have four poles for 

C 2 ± 2jC 
c 2 - 1 

However, before integrating in k we have to extend the range of integration from — oo to +00, 
which given that the whole integrand is even in k, is straightforward. Then, in sin(kr) we can 
keep just one of the two exponentials — for instance we can perform the replacement sin(A;r) — > 
sin(A;r) — icos(kr) = —ie tkr , which is allowed because of parity. We get 

G^V = 0) = fdfc fdC ^ • (71) 



e = . (70) 



7T^ r / j ' Cl> 



Now, since we are interested in studying the propagation of superluminal signals, we can restrict 
to the r > t case. That is, we study the retarded Green's function outside the usual relativistic 
light-cone. In this case, we can close the contour in the upper half-plane of k. There, at large k 
and fixed real ( the integrand behaves as 



e 



kr e ±i^C, 2 +k 2 t ^ ^2) 



which decays for r > t. Therefore we get no contributions from infinity, but we do get two pole 
contributions, from 



K = C\/^M = «(C) (73) 



l-2i/ C 
c 2 - 1 

which for real ( are both in the upper half-plane. Defining 



-C\/^r-^ = «(-C) (74) 



'c 2 + 2i/( 
c 2 - 1 

we have 



/3(C)^C\/-^-P, (75) 



G rct (r>t, j/ = 0) = 4(c2 _ i)7r2 - / rfC 



Ce^ S ^|P + (C--C) 



(76) 



Given the parity of the integrand in (, we can keep either term in brackets, modulo an overall 
factor of two. We will keep the first. As before, the factor sin(t (3(()) / j3(() is analytic in the whole 
(■-plane. On the other hand, «(C) has two branch-points, 

C = 0, and C = -2i. (77) 
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Our choice for the square root corresponds to making the branch-cut run between the two branch- 
points along the imaginary axis, as shown in fig. |6]a. 

If r > ct, we can close the contour in the upper half plane, where the integrand is analytic. 
The exponential behavior at infinity is 



e ira(0 e ±it(3(0 _ exp f £ . ( r ± ct )) 

which decays in the upper half plane as long as r > ct: 

G iet (r>ct,y = 0)=0 . 



(78 



(79) 



This shows that — as we expected — we can have superluminal signals traveling along the boundary, 
but they are confined to a new light-cone, with aperture determined by c. 

As a final check that the system is well-behaved, we must make sure that inside this new light- 
cone the retarded Green's function is not exponentially increasing in time — which would signal 
an instability of the vacuum upon local perturbations. This is simple to check in the interspace 



between the two light-cones, i.e. for t < r < ct. In this case we can still use (76). We now have to 
split the sine into two exponentials, 



sin(/%) 



2ip 



,m 



-i/3t 



) 



30) 



and analyze the two contributions separately. However, now for each individual contribution the 
analytic structure of f3(Q matters, and we have to keep track of that as well. /?(£) has a branch-cut 
running along the imaginary axis between and —2i/c 2 . For c 2 > 1, this is entirely contained 
in the branch-cut of a((). Now, for the term ~ e ira( -O e itm we can 

still close the contour in the 

upper half-plane. The behavior at infinity is 



ex P (vcfci^' ( r + ct )) > 



[a) 



Figure 6: The integration contours for the integral in C, as described in the text. The branch points are 
at £ = and C = —2i. 
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which vanishes for positive r and t. We thus get zero from this term. For the other term we have 
to close the contour in the lower half-plane. We still get a vanishing contribution from infinity 
as long as r < ct — as we are assuming — but upon deforming the contour as in fig. [6^, we are left 
with an integral along the cut: 



G iet (t < r < ct,y = 0) 



1 



oc 



4(c 2 
i 

dww A 
o 



-2i 



l) 71 " 2 T Jq 

exp 



3 ia(C)r-i/3(C)t 



m 



l/w) 



\J c 2 — l/w 



(82) 



33) 



where in the second step we performed the change of variable £ = —liw. The question now is 
whether this integral yields an exponential in time or r, with a positive real part in the exponent. 
We can split the integration region into two intervals, < w < 1/c 2 and 1/c 2 < w < 1. In 
the former, both square roots in the exponent are purely imaginary, and we get an oscillatory 
contribution. In the latter interval the square-root that multiplies r is still imaginary, but that 
multiplying t is real. However, our definition for the square root has a non-negative real part in the 
whole complex plane. Therefore this contribution is oscillatory in r and exponentially decaying in 
t. 

In summary, we have shown that the retarded Green's function for tt 

i) vanishes for \y\ > t, i.e. there is no superluminality in the direction orthogonal to the 
boundary; 

ii) vanishes for y = and r > ct, i.e. we have superluminality on the boundary, but we still 
have a causal cone with a finite aperture; 

in) is oscillatory inside this new causal cone — at least outside the usual relativistic light-cone — 
i.e. the system is stable against local perturbations. 



These are exactly the properties that one would have guessed from 7r's e.o.m. (46) 



3.1 The shape of the Green's function 

To conclude this section and generalize our results above to generic positions in the bulk, 
we wish to find the shape of the support of the retarded Green's function, with source centered 
at the origin. In the geometric optics approximation, its boundary at time t, as a function of x, 
will be the largest y value a light ray starting at the origin can reach given time t. We checked 
numerically that indeed the shape determined this way is correct. As usual, the causal structure 
is determined by the high-frequency limit. 

Let a light ray start at the origin and travel along the brane for a time t\ at the speed c, 
reaching the point x = dt\. It then leaves the brane and travels to point (x,y) along a straight 
path at a speed of 1 (see figure [7]). To find the optimal time t\ for the light ray to leave the brane, 
we wish to maximize y, given x and fixed total travel time t, as we vary over t\. We have 

y 2 + (x - chf = (t- ti) 2 , (84) 
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speed=c 



X 



Figure 7: Light ray traveling from the origin. 



which has a maximum at 



cx — t 

^1' 



ct — X 



ex— t 
c 2 -l> 

at 



Thus to reach the largest y value, the light ray should travel along the brane for a time t± 
and then leave the brane and make a bee-line for the point (x, 77=?^; )■ We see that t\ 
x = -, so at this value of x, it becomes most favorable to leave the brane immediately. 

The shape of the Green's function is shown in figure [8j From x = - c to x = ct, where it is most 



favorable to travel down the brane a bit before leaving, it is the straight line y 



ct—x 



From 



x = to x = -, where it is most favorable to leave the brane immediately, it is just the circle, 
y 2 + x 2 = t 2 . The circle and the line intersect at x = - and have the same slope there. This 
double teardrop shape represents the boundary of the support of the retarded Green's function. 
It expands linearly in t. It is of course rotationally symmetric among all of the x coordinates. 



4 Discussion 



Superluminality on a non-Lorentz invariant background does not imply per se a breakdown of 
causality in an effective field theory. It is nevertheless a troublesome signal which may indicate 
the absence of a relativistically invariant UV completion of the theory [12]. In this paper we 
showed that the DGP model exhibits superluminal propagation in otherwise physically sound 
backgrounds, that are exact solutions to the DGP equations. By computing small fluctuations of 
the scalar mode around exact solutions, we made sure that the pathology found here, which is 
of the same nature as that found in [121 02], is not an artifact of the decoupling limit M4 — > 00, 
A3 ~ M5/M4 = const. The scalar mode couples locally to physical sources, namely to the trace of 
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k y t 




Figure 8: Shape of the support of the Green's function at time t. 



the stress-energy tensor, but its equations of motion are non-local away from the strict decoupling 
limit. To be sure that the maximum speed of propagation we found using the local limit was a real 
effect, we computed the retarded Green's function of the scalar mode and found that on the 4D 
brane it is non-vanishing in an enlarged light cone. We finally showed that within this enlarged 
light cone, the Green's function is free of instabilities. This last check is important because it 
proves that the effect we found is truly a maximum speed of propagation of signals that exceeds 
c, instead of being a tachyon-like instability of the background. 

If we take the mildly conservative point of view that any theory that possesses a relativistically- 
invariant, causal completion should not allow superluminality, then we can use the results we found 
to rule out significant parts of the parameter space of DGP. Indeed, already the simple analysis 
presented in this paper shows that superluminality occurs whenever the 4D stress-energy tensor 
is a pure cosmological constant = —pg^ ui irrespective of the value of p. Of course our analysis 
relies on the existence of a domain of validity for the DGP model. DGP is after all an effective 
theory which gets corrected by higher-dimensional operators at length scales L5 ~ Mg -1 j^] To carry 
on our analysis we need of course Ldgp ~ Mf/Mf ^> L 5 . So, one way to escape the results of 
this paper is to have a 4D Planck mass M4 not too much larger than M5. This is what one finds 
by trying to construct DGP as the low-energy theory of a set of a few D-branes (or other kind of 
solitons) embedded in a higher-dimensional string theory. On the other hand, a large hierarchy 
M 4 ~ NM S 3> M 5 = M s /g1^ (with M s the string scale and g s the string coupling) might be 

5 Strictly speaking, there is a longer strong-coupling scale in the n sector of the theory, A^ 1 ~ M4/-M5 3> £5 
[6 . Nevertheless, classical solutions with (extrinsic) curvature in the strong-coupling regime are not necessarily 
outside the regime of validity of the effective theory [8] . 
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achieved by considering a theory with O(N) D-branes when N 3> l/gl^. Indeed, a hierarchy 
of scales between M\ and the UV cutoff in a theory with many species is a universal result that 
holds well beyond the DGP context [18] . In this particular realization the DGP scale is 

M\ N 2 g 2 

W ~ Wi - -jjf ■ (86) 

For our analysis to be valid we need L DGP 3> l/M s , that is 

Ng s > 1 . (87) 

In this case, the effective DGP action breaks down whenever curvatures are of order of the string 
scale. The tension of the D-brane arrangement — the four-dimensional cosmological constant en- 
ergy density, p, — has a classical contribution of order NMf/g s , and a quantum-mechanical one of 
order N 2 Mj. However for the latter we may hope for cancellations due to supersymmetry, so let 
us ignore it. The intrinsic curvature of the D-brane system can be small in our limit 

P Ml 



ft 



M 2 Ng s 



Still, the extrinsic curvature is always large, 

K »»~{q~ N 9s M s > M s , (89) 

thus signaling a breakdown of the effective field theory. In estimating the extrinsic curvature we 
used Israel's junction condition. In principle, there could be a cancellation between the localized 
DGP term and the cosmological constant, but in general we do not expect such a cancellation of 



the leading contribution (89). Indeed, such a cancellation does not take place for the solutions 
discussed in sects. 2.1, 2.2, where the extrinsic curvature is always of order p/M|. This is good, 
because those solutions are stable and have superluminal excitations — which we do not expect in 
a consistent D-brane arrangement in string-theory. 

In conclusion, in this simple construction there is no choice for the values of iV and g s yielding 
a consistent DGP effective theory, where both the DGP crossover scale and the curvature radii 
are longer than the cutoff. There is one dimensionless parameter, Ng s , controlling both how much 
longer than the cutoff L DGP is, and how much shorter than the cutoff the radius of curvature is. 
The results of this paper suggest that this obstruction may be generic and model-independent. In 
other words, it may be that the parameter space of DGP-like models that can be obtained from 
UV-complete theories is smaller than the space determined by demanding self-consistency of the 
effective field theory. Within that subspace, we expect to find no solution with pathologies like 
superluminal signal propagation. 
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